Level spectroscopy of the square-lattice three-state Potts model with a ferromagnetic 

next-nearest-neighbor coupling 



Hiromi Otsuka, Koutaro Mori, Yutaka Okabe, and Kiyohide Nomura^ 
Department of Physics, Tokyo Metropolitan University, Tokyo 192-0397 Japan 
^Department of Physics, Kyushu University, Fukuoka 812-8581 Japan 
(Dated: February 2, 2008) 

We study the square-lattice three-state Potts model with the ferromagnetic next-nearest-neighbor 
coupling at finite temperature. Using the level-spectroscopy method, we numerically analyze the 
excitation spectrum of the transfer matrices and precisely determine the global phase diagram. 
Then we find that, contrary to a previous result based on the finite-size scaling, the massless region 
continues up to the decoupling point with Z3 x Z3 criticality in the antiferromagnetic region. We 
also check the universal relations among excitation levels to provide the reliability of our result. 
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The concept of the Gaussian universality class pro- 
vides a paradigm for the unified descriptions and under- 
standing of the low-energy and long-distance properties 
of onc-dimcnsional (ID) quantum and two-dimensional 
(2D) classical systems. Its representation is given by 
the free-boson fixed point model, and its criticality by 
the conformal field theory (CFT) with the central charge 
c = 1. The antiferromagnetic three-state Potts (AF3SP) 
model on the square lattice is one of those to exhibit the 
Gaussian criticality at zero temperature 0, While 
the ground-state 0, B 0, 0, [S El an d finite-temperature 
properties were fully clarified, there still exists the 
considerable interest in its related extended-type mod- 
els 0, IS H- For instance, in 0, the emergence of the 
Z2 and Z3 criticalities were studied by introducing the 
AF3SP model with a staggered polarization field. 

In this paper we treat another model with a rather 
straightforward extension, i.e., the three-state Potts 
model with ferromagnetic next-nearest-ncighbor cou- 
pling defined on the square lattice A, whose reduced 
Hamiltonian 7i{K\, K2) = Rjk^T is given as 



U{K X ,K 2 ) = Ky S ai , n - \K 2 \ Y 5 aj , 

(j-k) b'.fe] 
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The first and the second sums run over all nearest- 
neighbor and ncxt-nearest-neighbor pairs, respectively, 
and the ternary variables dj take the values 0,1,2 (j G A). 
For the AF case (K\ > 0), this model is thought to be in 
the same universality class as the ferromagnetic six-state 
clock (F6SC) model Theoretical investigations in- 

cluding numerical ones have bee n p erformed to clarify 
the global phase diagram @,ELl2- However, its pre- 
cise estimation is not available. This is mainly because 
two types of Bcrezinskii-Kostcrlitz-Thouless (BKT) tran- 
sitions take part in the phase diagram, and the nu- 
merical methods used so far were insufficient to treat 
the BKT transitions. For instance, the phenomenolog- 
ical renormalization-group (PRG) method has been fre- 
quently used for the determination of second-order tran- 
sition points. However, as pointed out by several authors, 
it fails to estimate the BKT points ^S]- Therefore, a new 



strategy should be employed. In investigations of ID 
quantum systems, the level-spcctroscopy method which 
takes logarithmic corrections into account has been used 
for precise estimates of the BKT points, and its possible 
application to the 2D classical system was also discussed 
1141 . Therefore, we shall reconsider this long-standing 
model Q with a new methodological strategy in order 
to clarify the global phase diagram. While the detail is 
given in the following, newly obtained data agree with the 
exact results available in some limits of model parame- 
ters, and show that the phase diagram is constructed 
by the crossover of the criticality from the decoupling 
point (K 1 ,K 2 ) = (0,ln(l + \/3)) with Z 3 x Z 3 critical- 
ity (c = |), which is in contrast to the previous result 
based on the naive finite-size-scaling method 0. As a 
result, the precise determination of the phase diagram 
around the point turns out to be crucial to understand 
this model, whereas the Coulomb-gas picture provides a 
quantitative description only around the K\ = +00 line. 

Here we shall briefly refer to some relevant researches. 
First, the ground state of the AF3SP model {K\ = +00, 
A'2 = 0), which is equivalent to the ice-point six- vertex 
model and is thus critical 0, becomes off-critical at fi- 
nite temperature ^S|. It was pointed out that there are 
two types of excitations controlled by the thermal scaling 
field u = e 1 , i.e., the relevant and the marginal ones, 
and that both of these are necessary to explain the exotic 
corrections to scaling observed in the Monte Carlo data 
0. The energy operator ej = Y,/j,k) ^j,^k ( tnc scaling 
dimension x e = |) is the relevant excitation and brings 
about the correlation length £ oc u~ 1 ^ 2 ~ x ^ in its leading 
form. According to Delfino, the dual sine-Gordon La- 
grangian density in term of the bosonic field and its dual 
field <j>, 6 G [0,2tt/\/2] can describe this transition: 
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where y and y oc u are the coupling constants and a is 
a short-distance cutoff 0, ■ Since the Gaussian cou- 
pling K is \ for the AF3SP model at zero temperature, 
and thus the second term, the Zg symmetry-breaking 
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field, with the scaling dimension 6 is highly irrelevant 
to this transition, we can drop it for the discussions of 
the AF3SP model. However, it can become relevant for 
nonzero K 2 > (see below). 

Next, let us consider the case with u = (K\ = +00), 
where Eq. Q is equivalent to the exactly solved F model 
[T(| . Define v = e K2 ; then the v dependence of the Gaus- 
sian coupling is given as 1/9K = 1 — (I/71") cos - 1 (v 2 /2 — 1 ) 
[l7j . Thus the point v 2 = 2 (i.e., K = |) separates the 
massless (v < v 2 ) region and the massive (v > v 2 ) region 
with sixfold degeneracy, and the latter is stabilized by the 
relevant term y cos 6\^2(j). Further, the nonlinear terms 
in Eq. (J2J) are both irrelevant in the region V\ < v < v 2 
with vi = v/2 + 2cos(57r/9) ~ 1.2856 (K = J), so the 
critical phase spreads at least for small u and connects 
to the region in the limit in |. This argument is 
consistent with the so-called cell-spin analysis [l8|: Park 
showed that Eq. (|TJ with K 2 = Kx/2 is approximately 
reduced to the F6SC model in terms of the cell spins 
representing the unit cells of the sixfold ground states. 
Therefore, as a member of the F6SC universality class 
[T^ . the present model |JIJ is also expected to possess 
two BKT points vi(u) and v 2 (u), which is described by 
the dual sine-Gordon field theory Eq. 0. 

We shall first consider the system around v\ (u) where 
COS6V20 is irrelevant, and thus it is described by 

A = ^(W) 2 + ^c OS V20, K*\. (3) 

One of the authors (K.N.) investigated the properties of 
marginal operators on the BKT line: Especially, in this 
case, (1/ K) (V</>) 2 and \[2 cos y/W (the bosonized expres- 
sion of £j) hybridize along the renormalization-group flow 
and result in two orthogonalized operators, i.e., the U M.- 
like" and the "cos-like" operators 0|. Writing the for- 
mer and the latter as Oo and Ox, and defining the sys- 
tem on A with M (— ► 00) rows of L sites wrapped on a 
cylinder, near the multicritical point, their renormalized 
scaling dimensions are given as xq(1) ~ 2 — yo (l + |i) 
and X\{1) ~ 2 + yo (2 + where I = InL, (yo,yi) = 
(1/2K — 2,y), and the small deviation from the BKT 
point t = yi/yo — I. Since these operators are parts of 
Ci , the corresponding excitations possess the same sym- 
metry with the ground state (see below). Another im- 
portant operator is a relevant one: 2 = V2 sin 3y/2(j> 
whose dimension is expressed as x 2 (l) — (2 — yo) in 
the same region. Consequently, the level-crossing condi- 
tion xq(1) = ¥fX 2 (l) offers a finite-size estimate of the 
BKT point Vx(u,L) [3. 

Next we consider a region near v 2 (u) where cos \/29 is 
irrelevant. The effective Lagrangian is then given as 

£ 2 = ^(V0) 2 + ^ cos 6^20, *4 (4) 

In this case, the marginal operators (1/K) (V</>) 2 and 
V2 cos 6-^/20 also hybridize and result in the A'f-like and 
cos-like operators, Oq and 0\. Then, their scaling dimen- 
sions near the transition point are xq(1) ~ 2 — yo (l + |t) 



anda;i(Z) ~ 2 + y Q (2+ fi). Here, (y ,yx) = (18K-2,y) 
have been redefined. The operator 2 has the dimen- 
sion x 2 (l) ~ j [2 — yo (1 + 2t)\. Thus, the level-crossing 
condition for v 2 (u, L) can be given by Xq{1) = Ax 2 (l). Al- 
though there may exist some other level-crossing condi- 
tions for the determination of the BKT points, the levels 
focused on here and used in the following are relatively 
easy to be accessed by numerical calculations poj . 

Now, let us consider the system on A with M (— > 00) 
rows of L (an even number) sites, where j € A is specified 
by I G [1, L] and m € [1, 00]. For this system, we can de- 
fine the transfer matrix T(L) connecting (<7i. m , • • • , o\L,m) 
to (<7i, m +i, • • • , o"L,m+i), and denote its eigenvalues as 
{X n (L)} or their logarithms as {E n (L) = — In | A„(L)|} (n 
specifies a level). Then CFT provides direct expressions 
for c and x n (the scaling dimensions) of critical systems 
as E g (L) ~ Lf - ttvc/6L + b/L 3 and AE n {L) ~ ^x n , 
where E g is the smallest one corresponding to the ground 
state and AE n (L) = E n (L) — E g (L) is an excitation 
gap [U /, and b are the geometric factor 

(v = 1), the free energy per site, and a nonuniversal 
constant, respectively While we can calculate the renor- 
malized scaling dimension x n (l) from the excitation gap 
as AE n (L)/( 2 j^-), the symmetry properties are impor- 
tant for the specification of relevant excitations 0. In 
addition to the translation T and the space inversion "P, 
the S 3 symmetry associated with the global permutations 
of the ternary variables is also possessed by the present 
model JIJ. By defining the clock variable as Sj = e l27I " cr i/ 3 ; 
the two generators of the S3 group are expressed as the 
cyclic permutation S: Sj 1— > e l27r / 3 Sj for the Z3 sym- 
metry and the charge conjugation C: Sj 1— > s* for the 
charge conjugation symmetry (C). These operations are 
then related to the changes of the phase variable as S: 
V2(/> i-> V2(f> + 2tt/3 and C: V2<j) 1-* -V2(j) @. Appar- 
ently the marginal operators defined above are S3 invari- 
ant, and thus the corresponding excitation levels AE n (L) 
should be found in the subspace specified by this symme- 
try property. On the other hand, the relevant operator 
2 is Z3 invariant, but is odd for the charge conjugation. 
Further, since the sublattice symmetry property is deter- 
mined by the Z3 charge @, 2 is odd for this symmetry 
operation. Therefore, we calculate the excitation levels 
according to the indices for these symmetry operations. 

The exact diagonalization calculations of T(L) are car- 
ried out for systems with L =8-16. In Figs.^a) and[IJb), 
we plot examples of the K 2 dependencies of the scaled 
gaps x n = AE n {Kx,K 2 ,L)l(i^-) (or values multiplied 
by constants) at K\ = 1. Then, we find the points [i.e., 
the finite-size estimates vi^ 2 (u,L)] at which one of the 
above conditions is satisfied. Here, it is worthy of note 
that, in addition to the logarithmic correction, there is 
another type of correction between a lattice system and 
a continuous model which stems from the x = 4 irrele- 
vant operator L- 2 L- 2 1 in terms of CFT j2^. Therefore, 
we shall extrapolate them to the thermodynamic limit 
according to the least-squares fitting of the polynomi- 
als in 1/L 2 . In Table |U we give vi t2 (u,L) at K\ = 64 
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FIG. 1: Examples of the K2 dependence of scaling dimensions 
AE n (Ki,K 2 ,L)/(^-) at Ki = 1. (a) [(b)] shows x and 
16x2/9 (xo and 43:2). The correspondences between marks 
and system sizes are given in these panels. The crossing points 
give finite-size estimates v\,2{u,L). 



TABLE I: The finite-size estimates of the BKT points for 
systems close to the F model [1)1,2(0, L)]. Extrapolations to 
L — > 00 are performed according to the least-squares fitting 
of the polynomials in 1/L 2 . 



L 



10 



12 



14 



16 



«i(0, L) 1.31465 1.30449 1.29871 1.29517 1.29287 1.285 
v 2 (0,L) 2.02507 2.02626 2.02541 2.02404 2.02262 2.013 
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FIG. 2: The phase diagram of the model The open di- 
amonds (squares) with a curve exhibit the BKT line vi(u) 
[v2(u)]. The filled diamond and the square on the v axis cor- 
respond to (0, vi) and (0,1^), respectively. The filled circle at 
(1, l + \/3) denotes the decoupling point with Z3 x Z3 critical- 
ity. The open triangles with a curve exhibits the second-order 
transition line ^(m). The thick dotted line shows the first- 
order phase transition boundary. 



(i.e., a case almost equivalent to the F model) in order 
to compare our data with the exact limiting values vi_ 2 - 
While considerable size dependencies are visible in the 
estimates, we can check that the extrapolated data do 
not deviate more than 0.7% from the exact ones. 

In Fig. we give the phase diagram of our model Q 
in the space (u,v) — {e~ Kl , e K2 ). The open diamonds 
(squares) with a curve exhibit the phase boundary v\ (u) 
[i>2(u)], which separates the critical region from the disor- 
dered phase [from the ordered phase "Ordered (I)" with 
sixfold degeneracy]. The filled diamond and the square 
on the v axis show (0,^1) and (0, ^2), respectively. The 
filled circle at (u, v) — (1,1 + -^3) denotes the decou- 
pling point with Z3XZ3 criticality j24|. To complete the 
phase diagram, we also calculate the second-order phase 
transition point v% (u) with the ferromagnetic three-state 
Potts (F3SP) criticality in the ferromagnetic region u > 1 
(i.e., triangles) by employing the PRG method [3- The 
ordered phase "Ordered (II)" has threefold degeneracy, 
and the thick dotted line shows the first-order phase tran- 
sition boundary. From this figure we can find the follow- 
ing. While the cusplike behavior of the phase diagram 
is observed around the decoupling point [25j . two BKT 
lines vi^(u) do not merge into a single curve contrary to 
the observation in the previous work based on the PRG 
calculation Therefore, we conclude that the massless 
critical region continues up to the decoupling point and 
that neither a critical end point nor a first-order tran- 
sition line exists in the antifcrromagnctic region u < 1. 



Further, contrary to the previous calculation, the u de- 
pendence of v 2 (u) is monotonic and does not show a min- 
imum around u ~ 0.3. As a result, the phase diagram 
obtained by the level-spectroscopy method is consider- 
ably different from the previous one, and shows that the 
Gaussian critical region is stabilized in a wider area of 
the parameter space. 

Before moving to discussion, we shall check the consis- 
tencies among excitation levels on the BKT points. Since 
the amplitudes of logarithmic corrections are given by the 
operator- product-expansion coefficients, some universal 
relations among the scaling dimensions have been discov- 
ered 0H1: For instance, [2x (l) +x 1 (l)]/3 ~ 2 at vi(u) 
and [3x 2 (l) + x^(l)]/A ~ \ at v 2 (u). Here cc 3 (Z) is the di- 
mension of O3 = \/2 cos 3v / 2</> [the bosonized expression 
of the staggered polarization Pj = (-l) J E[ 3 ,fc]^jjJ 
which is S3 invariant and odd for the sublattice sym- 
metry HI23. In Table HU we give the scaling dimensions 
at K\ = 1. Although xo(l) and x\(l) [x 2 (l) and Xs(l)] 
considerably deviate from the value for the free-boson 
case 2 (^) due to the logarithmic corrections, their main 
parts cancel each other, so the average x av = (2xq+Xi)/3 
[^av = (3x2 + X3)/i] takes a value close to 2 (i). The 
consistency checks can be satisfied only if the systems 
are on the BKT lines, and the investigated levels based 
on symmetry properties have the expected physical inter- 
pretations. Therefore, Tabic UTI is helpful to demonstrate 
the reliability of our phase diagram. 

Finally, we discuss the nature of the phase diag ram. 
As the real-space RG calculations predicted [l2ll25l] . the 
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TABLE II: Examples of the L dependences of the scaling 
dimensions xq and X2 and the averages x av and x av (see text) 
on the BKT points «i,2(w) (Ki = 1). We extrapolate the 
finite-size estimates to L — > 00 using the least-squares fitting 
of the polynomial in 1/L 2 . 

L 8 10 12 14 16 00 

x (Z) 1.48366 1.52857 1.56357 1.59141 1.61404 — 

x av (l) 1.78633 1.83948 1.87883 1.90716 1.92759 2.006 

x 2 (l) 0.38137 0.38913 0.39454 0.39851 0.40155 — 

x av {l) 0.50271 0.50078 0.49933 0.49811 0.49704 0.492 



decoupling point (1, 1 + v^) with Z3XZ3 criticality (i.e., 
c = |) becomes unstable against relevant competing per- 
turbations and exhibits crossovers to the behaviors con- 
trolled by the critical fixed points with lower symmetries 
[28j . One of those perturbations is the energy opera- 
tor with scaling dimension | . For the ferromagnetic case 
u > 1, another one may be a product of order parameters 

on each sublattice which has the dimension 4z- Thcrc- 

15 

fore, the crossover exponent </>p = (2 — yg)/(2 — |) = ^ 
can predict the shape of the boundary around the point, 
and v 3 (u) agrees with this [3J. Further, we obtain, for 
instance, c ~ 0.800 at 1)3(11) (K\ = —0.25), so the bound- 
ary can be regarded as a massless RG flow to the F3SP 
fixed point with c = |. On the other hand, we also 
analyze vi^{u), and similar power law behaviors can be 



recognized. Although the estimated exponents take close 
values to ^ , the linearities of the log-log plots are rather 
lower than in the case of t>3(ii), and thus the estimates 
may acquire errors up to 10%. This may be due to the 
existence of the marginal operator in the intermediate 
critical region, which is absent in the ferromagnetic case. 
Since the central charge takes the values, e.g., c ~ 1.007 
and 1.009 at vi(u) and ^(u) {K\ = 1), respectively, the 
critical region can also be seen as a realization of the 
crossover to the Gaussian fixed point with c = 1 [2(| • 
However, its precise numerical analysis is still difficult, 
so this issue will be left as a future problem. 

In summary, we investigated the square-lattice three- 
state Potts model with the ferromagnetic next-nearest- 
neighbor coupling. The phase diagram obtained here re- 
veals new features that could not be found in previous 
research. With the use of the level-spectroscopy method, 
we can break new ground for studying 2D classical sys- 
tems with complex interactions. For more detailed study, 
we are now performing the Monte Carlo simulations, and 
we obtain preliminary data consistent with the present 
result of the phase diagram; they will be given in a forth- 
coming article. 
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